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ABSTRACT 


Effects  of  static  attitude  and  three  different  motions  on  vortex  behavior  and  its 
breakdown  over  delta  wings  have  been  investigated.  These  motions  include  1)  harmonic 
oscillation  or  ramp-and-hold  motion  in  roll  at  constant  resultant  angle  of  attack,  2) 
harmonic  oscillation  or  ramp-and-hold  in  pitch  and  3)  coning  at  constant  resultant  angle 
of  attack  A  basic  65°  delta  wing  configuration  was  tested  in  two  wind  tunnels  as  well  as 
in  a  water  tunnel.  A  roll  rig,  pitch  rig  and  coning  rig  (OPLEC)  were  employed  for  the 
investigation. 

Quite  different  flow  characteristics  have  been  found  within  the  spiral  vortex 
breakdown  region  even  under  steady  conditions.  The  point  where  vortex  breaks  down 
into  large  scale  turbulence  has  a  very  long  response  time  to  perturbations  and  is  the  major 
cause  for  the  dramatic  time  dependence  observed  in  the  airloads.  In  addition,  bi-stable 
vortex  breakdown  locations  are  observed  in  the  case  of  slender  delta  wings  when 
breakdown  is  located  close  to  the  trailing  edge.  Under  dynamic  conditions,  the  behavior 
of  the  vortex  and  its  breakdown  region  have  been  found  to  be  highly  dependent  on 
motion  waveform  and  reduced  angular  rates 

An  analytical  model  based  on  the  Nonlinear  Indicial  Response  Method  and  a 
formulation  suggested  by  Tobak  et  al,  and  capable  of  reflecting  the  above  effects  is 
proposed  to  predict  breakdown  location  under  dynamic  conditions. 

SYMBOLS 

b  wing  span 

c0  wing  root  chord 

D  vortex  core  diameter 

s  local  semi-span 

Uoc  freestream  velocity 

Xvb  chordwise  vortex  breakdown  location 

xvb  non-dimensional  vortex  breakdown  location  =  — ^ 

a  angle  of  attack 

a  non-dimensional  pitch  rate  =  -^A 

P  sideslip  angle 
8  deviation 
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A  standard  deviation 
A  leading-edge  sweepback  angle 
o  sting  angle  or  resultant  angle  of  attack 
(f>  aerodynamic  roll  angle 

<j>  non-dimensional  roll  rate 

x  convection  time 
co  angular  pitch  rate 
Q.  coning  rate 

Q.  non-dimensional  coning  rate 

INTRODUCTION 
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High-performance  combat  aircraft  should  be  capable  of  flying  at  high  incidence  and 
high  angular  rates,  where  flows  are  usually  dominated  by  strong  wing  or  LEX  leading- 
edge  vortices  needed  to  generate  large  aerodynamic  forces.  However,  the  behavior  of 
these  vortices  and,  especially,  their  breakdown  over  lifting  surfaces  are  the  major  cause 
for  the  observed  airload  nonlinearities  and  time  dependence  in  the  advanced  maneuvering 
regime.  Thus,  a  clear  understanding  of  the  behavior  of  vortices  and  their  breakdown  is 
an  essential  requirement  for  solving  the  flight  mechanics  problems  in  that  regime. 

To  date,  most  of  the  experimental  dynamic  studies  have  been  limited  to  motions  in 
one  plane,  either  pitch,  yaw  or  roll  .  Given  that  the  behavior  of  a  vortex  and  its 
breakdown  are  nonlinear  processes,  previous  linear  assumptions,  that  allow  the 
superposition  of  the  effects  due  to  the  motion  orthogonal  components  is  not  legitimate. 

In  the  nonlinear  regime,  where  vortex  behavior  and  its  breakdown  depend  on  all  the 
motion  variables  and  their  histories,  the  response  to  motions  in  one  plane  clearly  will  be 
influenced  by  the  presence  of  other  motions.  The  practical  difficulties  of  imparting  all 
motions  at  once,  induced  Tobak  and  his  colleagues  to  investigate  this  coupling  in  a 
nonlinear  context1.  They  found  that  for  nearly  rectilinear  flight  paths  in  the  non-linear 
regime,  the  aerodynamic  response  to  an  arbitrary  motion  may  be  obtained  from  the 
contributions  from  four  simple  motions  which  are  all  at  a  constant  inclination  of  the  body 
axes  with  respect  to  the  velocity  vector.  These  four  simple  motions  are:  1)  fixed  attitude, 
2)  oscillation  in  roll  at  constant  resultant  angle  of  attack,  3)  oscillation  in  pitch  about  a 
fixed  resultant  angle  of  attack  and  4)  coning  at  constant  resultant  angle  of  attack  (Fig.  1). 

Since  the  aerodynamic  coefficients  are  functions  of  the  flow  field,  which  at  high 
incidence  is  overwhelmingly  vortex  dominated,  it  stands  to  reason  that  a  similar 
methodology  as  used  by  Tobak  can  be  used  as  a  first  step  in  describing  vortex  behavior  in 
the  case  of  complex  motions.  Thus,  the  behavior  of  the  vortex  and  its  breakdown  under 
those  four  simple  motions  are  discussed  in  this  paper. 

EXPERIMENTAL  SET-UP 


Studies  of  the  effects  of  the  above  motions  were  conducted  at  the  IAR  6x9  ft,  USAF 
7x10ft  wind  tunnels  and  IAR  15x20  water  tunnel  respectively.  A  roll  rig,  a  pitch  rig23 
and  an  orbital-platform-coning  (OPFEC)  rig45  6  were  employed  in  the  studies  (Fig.  2  to 
Fig.  4).  The  roll  rig  is  capable  of  imparting  an  harmonic  motion  with  amplitude  and 
frequency  of  up  to  40°  and  18  Hz  respectively  with  a  maximum  roll  angle  offset  ±50°. 
The  pitch  rig  rotated  the  model  around  an  axis  5.39  in.  behind  the  trailing  edge  over  a 
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range  of  angle  of  attack  from  -1°  to  90°  with  minimum  transition  time  of  75  ms.  In  the 
case  of  the  OPLEC  apparatus,  the  model  is  mounted  on  a  sting  attached  to  an  annular 
platform  riding  on  the  outer  surface  of  a  stationary  cylindrical  test  section  insert.  This 
configuration  results  in  minimal  aerodynamic  interference  as  the  strong  support 
interference,  present  in  conventional  rigs  caused  by  the  rotating  support  arm,  is 
eliminated. 

The  basic  configuration  in  the  investigation  is  a  65°  delta  wing  (Fig.  5).  Static 
sweeps  as  well  as  large-amplitude,  high-rate  harmonic  and  ramp  and  hold  motions  in 
pitch  and  roll,  as  well  as  coning  motion  were  performed.  The  tests  were  conducted  at  a 
Reynolds  number  of  36000  and  3.6  million  per  foot  in  the  water  tunnel  and  wind  tunnel 
respectively.  In  addition,  some  slender  delta  wings  were  tested  statically  in  IAR  water 
tunnel. 

In  the  wind  tunnel  experiments,  leading-edge  vortex  was  visualized  by  seeding  the 
flow  with  smoke  or  natural  condensation  and  applying  a  laser  light  sheet  oriented  normal 
to  the  body  axis  thus  yielding  a  section  view  of  the  vortex  or  longitudinally  passing 
through  the  vortex  core  over  most  of  its  length.  In  the  smoke  tests,  the  plume  of  smoke 
was  made  to  impinge  near  the  nodal  point  of  attachment  on  the  windward  side.  If  the 
lighting  set-up  and  gray  scale  sensitivity  is  fixed  during  the  experiment,  the  gray  level  is 
approximately  proportion  to  w/ux.  Thus,  the  image  reflects  the  behavior  of  the  rotational 
vortex  core7.  In  order  to  obtain  a  reasonably  high  temporal  resolution,  images  were 
captured  by  a  high-speed  video  camera  with  a  frame  rate  of  up  to  1000  fps.  Video  frames 
were  individually  digitized  using  an  8  bit  gray  scale  (black=0,  white=255).  In  general,  an 
intact  vortex  appears  as  a  ring  of  high  brightness  with  a  dark  core  whereas  a  burst  vortex 
appears  as  an  evenly  illuminated  disk  with  somewhat  more  diffuse  edges.  Occasionally, 
the  darker  sub-core  is  not  discernible  due  to  saturation  of  the  video  camera.  The  edge  of 
the  vortex  core  was  assumed  to  be  where  the  maximum  gradient  occurs.  Once  the  edge 
was  defined,  its  diameter  and  position  relative  to  the  model  were  determined.  After  the 
center  of  the  vortex  was  obtained,  the  spanwise  gray  level  profiles  along  the  vortex  axis 
were  further  obtained.  The  criteria  for  determining  vortex  breakdown  in  the  section  view 
is  a  little  downstream  of  that  based  on  the  kink  point  of  the  vortex  filament  in  a 
longitudinal  view,  as  it  is  difficult  to  detect  the  kink  point  in  the  former  view.  In  the 
water  tunnel  experiments  dye  is  injected  through  two  0.008  in.  diameter  ports  located  at 
5%  c0  and  half  local  semi-span  on  the  windward  side.  Two  video  cameras  provide  a  top 
and  side  view  of  the  vortices.  Breakdown  was  defined  at  the  kink  of  the  vortex  filament. 

EXPERIMENT  RESULTS 
Basic  behavior  of  vortex  in  steady  situation 

As  described  by  Lambourne8  a  vortex  exhibits  three  stages  1).  Axial  filament  is 
almost  straight  up  to  a  kink  point.  2)  Downstream  of  this  point  the  vortex  remains  intact 
but  it  is  deformed  into  a  spiral  that  rotates  in  the  opposite  direction  as  the  spiral 
streamline  of  the  intact  vortex.  3)  After  a  few  turns,  the  vortex  eventually  breaks  down 
into  large-scale  turbulence.  It  is  worthwhile  to  mention  that  in  Lambourne’ s  definition, 
the  vortex  is  still  intact  downstream  of  the  kink  point  which  is  different  from  the 
prevalent  definition. 

These  three-stages  were  clearly  seen  in  the  experiments.  Large-scale  fluctuations 
were  found  in  the  breakdown  region  even  under  static  conditions.  As  examples,  Fig.  6 
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shows  two  typical  patterns  obtained  for  the  same  model  and  attitude.  To  facilitate 
descriptions,  the  vortex  filament  kink  is  denoted  as  point  A  in  Fig.  5  while  the  subsequent 
turbulent  breakdown  is  denoted  as  point  B.  Two  white  spots  marked  on  the  wing  surface 
are  employed  as  reference  for  measuring  distances.  The  distance  between  the  two  marks 
is  0.34  root  chord,  c0.  Interestingly,  the  breakdown  process  can  cover  almost  0.3  c0  as 
seen  from  Fig.  6.  More  importantly,  the  behavior  of  the  two  points  are  quite  distinct. 
Whereas  the  locations  of  point  B  are  nearly  the  same  in  both  cases,  those  of  point  A  are 
radically  different. 

The  high  level  of  unsteadiness  in  the  breakdown  region  is  illustrated  in  Fig.  7  where 
selected  pairs  of  consecutive  frames  are  shown.  The  frames  in  the  first  and  second  pairs 
are  separated  by  4  ms  while  the  frames  in  the  last  pair  are  taken  1  ms  apart, 
corresponding  to  0.60  and  0.15  convection  times,  c/UM,  respectively.  These  typical 
figures  show  that  the  length  of  the  spiral  breakdown  region  can  change  significantly 
within  a  convection  time  indicating  that  the  propagation  speed  of  point  A  can  approach 
0.5Ucos«  whereas  the  very  small  displacement  of  point  B  reflects  its  much  lower 

propagation  speed.  Furthermore,  Fig.  7  seems  to  suggest  that  for  a  given  angle  of  attack 
when  the  instantaneous  flow  condition  promotes  vortex  breakdown,  point  A  moves 
rapidly  forward  while  the  response  of  point  B  is  extremely  sluggish,  resulting  in  an 
increase  of  the  helix  angle,  a  reduction  in  the  spiral  diameter  with  the  pitch  remaining 
more  or  less  constant.  On  the  other  hand  when  the  instantaneous  flow  condition  delays 
breakdown,  point  A  moves  rapidly  downstream  with  point  B  again  responding  very 
slowly,  resulting  in  a  decrease  of  the  helix  angle  and  an  increase  in  the  spiral  diameter, 
with  pitch  again  remaining  more  or  less  constant.  The  length  of  the  breakdown  region 
therefore  changes  mainly  as  a  result  of  changes  in  the  number  of  spiral  turns  contained 
therein.  The  additional  loops  needed  to  lengthen  the  breakdown  region  seem  to  be 
generated  by  curling  up  of  the  previously  intact  vortex  just  forward  of  point  A.  Likewise 
a  straightening  of  the  upstream  turns  within  the  breakdown  region  seems  to  cause  its 
shortening. 

Energy  in  the  spiral  breakdown  region  decays  gradually  between  points  A  and  B 
resulting  in  a  slower  streamwise  reduction  in  suction  than  in  the  case  of  bubble 
breakdown.  Moreover,  the  much  slower  response  of  point  B  clearly  dictates  the  very 
important  time  lags  observed  in  the  breakdown  motion  reflected  by  the  movement  of 
Point  suggesting  that  Point  B  is  a  more  important  characteristic  parameter  than  Point  A. 

Furthermore,  it  has  been  found  that  the  breakdown  unsteady  behavior  is  significantly 
affected  by  geometric  discontinuities  on  the  model  surface  possibly  through  their  effect 
on  the  local  streamwise  pressure  gradient.  Significantly  larger  fluctuations  of  the 
breakdown  location  were  observed  when  breakdown  occurred  over  the  leading-  and 
trailing-edge  bevels  as  demonstrated  by  the  standard  and  maximum  deviations  plotted  in 
Fig. 8.  This  result  was  obtained  in  a  water  tunnel  for  a  65°  delta  wing  with  a  10°  leading- 
and  trailing  edge  bevels.  In  order  to  assess  the  stability  of  the  breakdown  process, 
experiments  were  conducted  in  which  an  obstruction  was  momentarily  introduced  as  a 
pulse  input  at  a  location  aft  of  breakdown  for  a  75°  delta  wing  at  a=36°,  34°  and  31.95° 
corresponding  to  xVB/c=  0.51,  0.67  and  0.94  respectively.  In  the  first  two  cases,  the  pulse 
input  resulted  in  a  momentary  forward  displacement  of  the  breakdown  location  which 
quickly  returned  to  its  original  position  once  the  obstruction  was  removed  as  depicted  in 
Fig.  9  and  Fig.  10.  The  breakdown  location  on  the  other  wing  half  exhibited  the  opposite 
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behavior  but  also  returned  to  its  original  position  after  removing  the  disturbance.  Similar 
responses  were  obtained  for  wings  with  a  smaller  sweepback  angle. 

In  the  case  where  breakdown  was  in  the  vicinity  of  the  trailing  edge  at  a=31.95°, 
however,  a  totally  different  phenomenon  was  observed.  Breakdown  initially  occurred 
symmetrically  at  xVB  -0.94;  and  if  the  obstruction  was  applied  sufficiently  aft  of 
breakdown,  the  response  was  analogous  to  those  at  higher  angles  of  attack.  However, 
after  introducing  and  withdrawing  the  obstruction  closer  to  the  breakdown  location,  the 
breakdown  location  unexpectedly  moved  forward  to  a  stable  position  5%  upstream  of  the 
initial  symmetric  location.  In  the  meantime  the  breakdown  point  of  the  other  vortex 
commenced  a  gradual  motion  that  would  place  it  well  aft  of  the  trailing  edge,  taking  over 
hundred  free-stream  convection  times  to  reach  a  stable  position  at  xVB~1.4.  Application 
of  the  disturbance  aft  of  the  latter  could  either  lead  to  a  reversal  of  the  situation  or  a 
symmetric  arrangement.  An  edited  time  history  of  the  breakdown  locations,  where  long 
waiting  times  have  been  removed  for  clarity,  is  shown  in  Fig.  11. 

The  various  breakdown  responses  to  the  above  phenomena  can  be  readily  discerned, 
which  means  that  more  than  one  stable  location  exists  for  the  time-averaged  breakdown 
point.  The  attractor  associated  with  the  symmetric  arrangement  was  smaller  than  that  for 
the  asymmetric  ones  if  the  breakdown  location  is  close  to  the  trailing  edge.  It  is, 
however,  worthwhile  to  mention  that  this  phenomenon  is  only  observed  over  more 
slender  wing  with  breakdown  location  near  the  trailing  edge. 

Motion  effect  on  vortex  behavior  and  its  breakdown 

Under  dynamic  conditions,  the  behavior  of  the  leading-edge  vortex  is  considerably 
more  complex  due  to  the  motion  effects  which  are  superimposed  to  the  steady  ones.  Fig. 
12  shows  the  non-dimensional  vortex  diameter  vs.  roll  angle  over  the  right  wing  half  at 
several  stations  for  different  rolling  oscillation  frequencies.  These  results  were  obtained 
from  gray  level  analysis  of  the  vortex  section  view.  The  static  results,  based  on  the 
average  value  over  the  intact  vortex  are  also  included  for  comparison.  Each  dynamic 
trajectory  has  two  branches,  the  left  and  right  ones  corresponding  to  the  upstroke  and 
downstroke  respectively.  The  first  point  on  the  left  branch  and  the  last  one  on  the  right 
branch  correspond  to  the  formation  and  breakdown  of  the  vortex  at  the  indicated  stations. 
For  a  given  roll  angle,  the  dynamic  D/s  is  larger  than  the  corresponding  static  value 
during  the  upstroke,  while  it  is  smaller  during  the  downstroke.  The  difference  between 
the  dynamic  and  static  values  increases  with  amplitude  and  frequency,  although  at  higher 
angular  rates  this  process  tends  to  saturate. 

There  is  abundant  evidence  that  the  motion  of  leading-edge  vortex  breakdown  plays 
an  important  role  in  creating  the  airload  non-linearities  and  time  dependence  observed 
under  dynamic  conditions.  Considerable  effort  has,  therefore,  been  directed  at  observing 
the  movement  of  breakdown  location,  which  must  be  done  by  means  of  an  averaging 
process  to  eliminate  the  above  fluctuations  a  process  which  is  greatly  simplified  by  using 
point  B  because  of  its  inherent  slow  response 

The  effect  of  rolling  oscillation  frequency  on  breakdown  motion  over  the  right  wing 
half  is  shown  in  Fig.  13.  The  response  to  the  rolling  oscillation  exhibits  a  hysteretic 
behavior  even  at  very  low  reduced  frequencies  and  becomes  more  pronounced  at  higher 
frequencies.  The  differences  between  the  results  at  k=0.14  and  k=0.2  are  relatively  small 
due  to  the  saturation  process  alluded  to  above.  The  longitudinal  range  traversed  by  the 
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breakdown  point  decreases  at  higher  frequencies  due  to  its  limited  propagation  speed, 
which  has  been  found  to  be  approximately  0.1  Uoo  and  0.2  Uoo  in  the  forward  and  aft 
directions  respectively.  Sections  of  the  loops  are  not  shown  due  to  insufficient  frames  to 
reliably  establish  the  breakdown  location.  The  effect  of  amplitude  is  shown  in  Fig.  14. 
In  this  case  the  loop  tends  to  expand  rather  evenly  with  amplitude  without  exhibiting 
saturation,  suggesting  a  more  linear  effect  than  that  of  frequency. 

Under  transient  ramp-and-hold  rolling  conditions,  significant  time  lags  are  similarly 
present  with  respect  to  their  corresponding  time-averaged  static  locations.  As  examples, 
the  instantaneous  vortex  breakdown  location  measured  during  ramp-and-hold  from  +16° 
to  -16°  roll  angle  or  reversed  from  -16°  to  +16°  motion  with  different  ramp  rates  are 
depicted  in  Fig.  15.  The  breakdown  location  shifts  more  as  the  ramp  rate  increases.  For 
example,  the  average  breakdown  location  was  delayed  by  1°  roll  angle  at  the  rate  of 
±0.00072,  while  it  was  delayed  4°  at  the  roll  rate  of  ±0.0028. 

The  effect  of  pitch  motion  on  vortex  behavior  can  be  observed  from  Fig.  16  to  Fig. 
18.  As  an  illustration,  the  time-averaged  position  for  the  case  of  pitch-down  transition 
from  a=30°  to  20°  are  depicted  in  Fig.  16a  where  the  static  ones  are  superimposed  for 
comparison.  Initially,  breakdown  is  at  its  static  location  and  moves  aft  when  the  angle  of 
attack  decreases.  However,  in  spite  of  the  small  non-dimensional  angular  rate,  the 
instantaneous  breakdown  location  rapidly  falls  behind  its  static  counterpart.  Examples  of 
the  effect  of  pitch-up  motion  between  a=20°  and  50°  at  different  angular  rates  are  shown 
in  Fig.  17a  and  Fig.  17b.  In  both  cases  breakdown  is  located  aft  of  the  trailing  edge  at  the 
motion  onset  and  move  on  to  the  wing  at  a  larger  angle  of  attack  dependent  on  the 
angular  rate.  Interestingly,  at  the  lower  rate  the  instantaneous  breakdown  location  tends 
to  approach  its  static  value  (Fig.  17a),  while  the  angular  rate  is  still  increasing.  Results 
corresponding  to  the  same  motion  but  at  a  higher  angular  rate,  show  that  breakdown 
reaches  the  trailing  edge  at  an  even  higher  angle  of  attack  (a=31°)  far  later  than  that 
under  static  conditions. 

In  addition,  pitch  motion  also  causes  axial  compression  and  expansion  of  the 
breakdown  region  between  points  A  and  B.  As  seen  in  Fig.  16  point  A  tends  to  “catch 
up”  point  B  due  to  its  higher  speed  of  response  during  pitch-down  while  the  distance 
seems  to  grow  reaching  a  maximum  value  at  a -31°  during  pitch-up  (Fig.  17).  This 
breakdown  behavior,  together  with  the  previously  described  unsteady  one,  appears  to 
suggest  that  although  point  A  tends  to  react  much  faster  than  B,  the  latter  in  some  way 
restricts  the  motion  of  the  former. 

Fluctuations  caused  by  flow  unsteadiness  are  also  observed  under  dynamic 
conditions.  Fig.  18,  as  an  example,  shows  two  consecutive  images  taken  at  instantaneous 
angles  of  attack  a=24.43°  and  a=24.40°  during  a  pitch  down  motion.  The  time 
difference  is  1  ms.  Point  A  in  Fig.  18b  should  move  aft  compared  with  that  in  Fig.  17, 
instead,  it  moves  significantly  forward.  In  the  meantime  there  is  no  significant  change  in 
point  B  consistent  with  the  previously  described  sluggishness  of  the  point  B  response. 

Coning  motion  effects  on  the  leading-edge  vortices  and  their  breakdown  are 
illustrated  in  Fig.  19.  As  examples,  the  measured  breakdown  locations  for  different 
coning  rates  and  resultant  angles  of  attack  are  depicted  in  Fig.  20.  Generally  speaking, 
the  coning  motion  promotes  breakdown  on  the  advancing  side  of  the  wing  while  delaying 
it  on  the  receding  side  (based  on  the  part  of  the  model  forward  of  the  intersection  with 
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coning  axis).  The  coning  motion  also  affects  the  trajectory  of  the  vortex  core  by  causing 
it  to  acquire  a  gentle  spiral  shape.  A  positive  coning  rate  (clockwise  in  back  view) 
displaces  the  trajectory  of  the  port  vortex  core  outboard  on  the  forward  portion  and 
inboard  in  the  aft  portion  of  the  wing  half. 

GENERALIZED  MODELING  OF  MOTIONS’  EFFECTS 

So  far  vortex  behavior  under  four  different  motions  have  been  studied  on  models 
with  the  same  geometry.  In  order  to  predict  vortex  behavior  due  to  an  arbitrary  motion  in 
the  nonlinear  regime  a  generalized  model  based  on  Tobak's  Nonlinear  Indicial  Response 
method  is  proposed.  Such  an  approach  becomes  necessary  because  in  this  regime  there 
are  couplings  between  the  effects  of  the  motion  orthogonal  components  and  moreover, 
these  effects  are  also  dependent  on  the  motion  history.  Obviously  the  linear  model  which 
assumes  that  the  indicial  responses  are  independent  of  motion  history  and  allows  the 
superposition  of  independent  effects  due  to  the  orthogonal  motion  components  is  not 
applicable  under  these  conditions. 

Thus  Tobak's  approach,  originally  developed  to  predict  airloads  is  adopted  and 
extended  to  a  general  variable  representing  the  flow  field  or  more  specifically  the  vortex 
breakdown  location-  at  high  angle  of  attack  and  high  angular  rate.  This  appears  to  be 
legitimate  as  aerodynamic  coefficients  reflect  the  integrated  effect  of  the  flow  field  which 
at  high  incidence  is  clearly  vortex  dominated.  Moreover,  experimental  results  have 
demonstrated  that  the  vortex  turbulent  breakdown  location  (point  B)  is  less  affected  by 
flow  unsteadiness  and  is  thus  considered  to  be  the  variable  of  interest,  which  as  a  first 
approximation  may  be  taken  as  a  continuous  function  of  the  motion  variables  except  in 
the  cases  discussed  above  of  highly  swept  delta  wings  when  breakdown  is  near  the 
trailing  edge,  or  the  model  leeward  surface  exhibits  geometric  discontinuities.  Given  the 
dominant  effect  of  breakdown  on  aerodynamic  loads,  if  its  movement  is  solved  it  can  be 
taken  as  an  internal  variable  towards  the  final  estimation  of  aerodynamic  loads  in  the 
nonlinear  regime.  Preliminary  studies  conducted  in  1992  and  later  in  2000  show 
prospective  results  which  agree  with  the  experimental  results910.  A  similar  idea  but  using 
a  different  approach  was  proposed  by  Goman  and  Khrabrov11 

According  to  this  method,  the  aerodynamic  response  (breakdown  location  in  this 
case)  can  be  built  up  from  the  responses  to  a  limited  number  of  well-defined 
characteristic  motions.  In  general,  xVB  can  be  expressed  as  a  functional  of  a(q),  4>f£,)  and 
fl(q),  in  an  arbitrary  motion  over  the  interval  -°«<^<t: 

xVB(t)  =  xVB[o(^),(|)(^),i2(^)]  (la) 

or  explicitly  in  the  form  of  an  indicial  response  functional 

xVB(t)  =  xVB[a(^),<K^),^©;t,T]  (lb) 

which  means  that  xVB(t)  is  both  a  functional  of  o(c,),  (xqj.tKc,)  and  a  function  of  t  and  x.  If 
the  functional  is  Frechet  differentiable  over  ( a(q),<^(q),  Q(c,)j  for  the  motion  at  hand,  the 
following  expression  applies: 

f'(o,  (]),£!)  =  fa(a,(j),  iDAo  +  f^io,  (]),  H)A(|)  +  f^(o,  <|),  tI)A£2  (2) 

where  /  stands  for  xVB  in  this  case  and  Ac ,  A0 ,  AQ  are  unit  step  changes  in 
respectively. 
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Based  on  the  NIR  methods  (Eqn  (19)  in  Ref.  1),  the  vortex  breakdown  location  can 
then  be  given  by 

X  VB  (t)  =  X  VB0  (0)  +  J  X  VB1  [a®,  <H£),  t,  x]  dx  +  J  X  VB2  [o(qj,  <K^),  H©;  t,  x]  dx 

o  o 

+  J  X  VB3  [<*©,  <K5).  t,  x]  dx  (3) 

0 

where  vVB1.,  xVB2  ,  vVB3  are  the  indicial  responses  to  step  changes  in  o,  (]),  Q.  respectively. 
Contrary  to  the  linear  case,  here  the  response  to  each  step  depends  on  the  motion  history 
prior  to  the  step  and  therefore  is  a  nonlinear  function.  This  superposition  is  valid  even 
though  the  system  is  nonlinear  if  the  equilibrium  flow  at  each  step  onset  is  asymptotically 
stable  to  small  perturbations. 

These  response  functions  could  be  further  understood  as  nonlinear  responses  to 
instantaneous  angular  rate  in  a  certain  period  of  motion  history.  Heuristically,  the  effects 
of  angular  rates  can  be  visualized  by  a  continuously  deforming  wing  geometry  whose 
surface  slopes  match  the  effective  local  “angles  of  attack”,  “sweepback  angle”,  etc., 
induced  by  the  angular  rates. 

Thus  the  modeling  problem  can  be  reduced  to  two  steps:  firstly,  the  determination  of 
the  equivalent  geometric  changes  due  to  angular  rates  in  order  to  establish  the  quasi¬ 
steady  breakdown  response  and  secondly,  the  determination  of  the  unsteady  response  to 
these  changes  including  motion  history  effects.  The  first  step  is  relatively  simple  and 
helps  to  understand  some  of  the  experimental  results.  The  second  one  is  far  more 
complex  and  has  so  far  been  successfully  completed  for  the  rolling  motion  case. 

In  the  pitching  case  the  equivalent  normal  surface  coordinates  can  be  found  from  the 
normal  velocity  distribution  to  be  (Fig.  21a). 


z(t) 


_  1  «(t)  2 

“  2  U„ 


(4) 


In  the  case  of  rolling  motion,  the  surface  normal  coordinates  are  obtained  by 
simulating  the  roll-rate-induced  chordwise  angle  and  section  angle  normal  to  the  leading- 
edge,  and  for  zero  roll  angle  are  given  by  8: 

z(t)=  [(— cor2A(t)  + - - - )|y|-x]  (5) 

coso(t)  sin2A(t)(l±tano(t)sin(])(t)tanA(t) 

where  ±  sign  correspond  to  the  right  and  left  wing  halves  respectively.  It  is  interesting  to 
note  that  for  (])=0  the  wing  surface  is  anti- symmetrical  about  the  body  axis  of  x  and  is 
characterized  by  straight  leading  edges  as  illustrated  in  Fig.  21b. 

Finally,  in  then  coning  motion,  the  lateral  velocity  component,  U  ,  induced  by  the 
coning  can  be  estimated  by  means  of  the  following  kinematic  equation: 

Uy  (x)  =  U  00sincj(t)  sin  (f)(t)  +  Q(t)xsino(t)cos(p(t)  =  U00sino(t)(sin(p(t)  +  2Q(t)(x  /b)coscp(t)) 

where  the  chordwise  coordinate  x  has  its  origin  at  the  coning  center.  This  additional 
velocity  component  results  in  changes  in  effective  sideslip  angle  and  additional  leading- 
edge  sweep  angle: 

(3(x)  =  sin_1[Uy  (x)/U«,]  =  sin-1[sina(t)(sin<|)(t)  +  2£2(t)(x/b)  cos  cf>(  t))] 


(6a) 
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AA(x)  =  tan  1(Uy/Ux)  =  tan  1[tano(t)(sin(|)(t)  +  2i2(t)(x/b)cos(|)(t))]  (6b) 

For  positive  coning,  the  lateral  velocity  component  and  sideslip  angle  is  negative  in 
the  forward  portion  and  positive  in  the  aft  portion  as  seen  in  Fig.  21c. 

Moreover,  the  local  aerodynamic  roll  angle,  cf>a,  is  also  affected  as  follows: 

(^a (x)  =  tan_1(Uy  /Uz) 

_  _  (6c) 

=  tan-1  { [tan  (|)(t)  +  2i2(  t)(x  /  b)]  /[I  -  2i2(t)(x  /  b)  tan  (|)(t)] } 

where 

Uz  =11^  sin cy(t)(cos (f)(t)  —  2il(t)(x/ b) sin <()(t)) 

From  a  Lagrangian  viewpoint  the  wing  has  different  camber  on  each  side  as  shown  in 
Fig.  2 Id.  At  (|)(t)=0  the  normal  coordinates  of  the  equivalent  surface  is  very  simple  and  is 
given  by 

z  =  20.— 
b 

All  those  changes  in  effective  leading-edge,  camber,  especially  in  the  forward 
portion,  affect  the  movement  of  the  vortex  breakdown. 

Other  than  the  breakdown  location,  coning  motion  also  affects  the  trajectory  of  the 
vortex  core.  When  the  wing  is  coning,  the  vortex  core  exhibits  a  spiral  pattern  as  seen  in 
Fig.  22b  and  22c  .  Since  coning  introduces  a  lateral  velocity  component  as  seen  in  Fig. 
21c,  this  lateral  velocity  component  will  push  or  pull  the  vortex  core  spanwise  according 
to  the  direction  of  the  lateral  velocity.  Thus,  for  example,  in  the  positive  coning  rate  the 
trajectory  of  vortex  core  moves  outboard  on  the  forward  potion  of  receding  wing  half 
while  it  moves  inboard  on  the  aft.  This  movement  will  largely  affect  the  rolling  moment. 

Given  the  extreme  complexity  of  obtaining  the  previously  mentioned  unsteady  breakdown 
response,  experiments  utilizing  models  deformed  according  to  the  above  criteria  can  be 
conducted  in  order  to  facilitate  the  process. 

As  an  example,  this  approach  has  been  successfully  used  in  the  case  of  rolling  motion 
yielding  the  following  expression.8: 

xVB(t)  =  Xs«j>(t))  +  Jtl_T  X  d  (t  —  x)<j>(  x)dx  (7) 

The  first  term  represents  the  effect  of  steady-state  at  time  t  while  the  second  one  is 
dynamic  effect  including  quasi-steady  and  unsteady  effect.  The  later  can  be  explicitly 
expended  as  the  summation  of  quasi-static  and  motion  history  effects. 

X VB  (t)  =  Xs  «Kt))  +  xq  «j>(t))  Xs«Kt))  +  Jt‘_T  Xu  (t  -  x)<j>(x)dx  (8) 

Through  quasi-steady  experiment12  with  deformed  wing,  roll  rate  induced  quasi¬ 
steady  effect  is  found  as  shown  in  Fig.  23.  The  corresponding  expression  is: 

Xq((j>(t),  a)  =  kqcj>(t)  =  tanttCO  ^>(t)  (9) 

From  experimental  results,  a  simplified  response  function  is  obtained: 

X  „(♦;«<*),  n<t»  = sin  ^ 


(10) 
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Fig.  24  shows  the  good  comparisons  between  measured  and  predicted  breakdown 
locations.  Similar  procedures  can  in  principle  be  applied  to  pitch  and  coning  motions. 
When  all  those  unknown  in  Eqn.  (3)  are  solved,  the  movements  of  vortex  breakdown 
location  under  arbitrary  motion  could  then  be  obtained. 
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STEADY  ANGLE  OF  ATTACK  PITCH  OSCILLATIONS 


Fig.  1  Characteristic  motions  in  the  aerodynamic  axis  system 


Fig.  2  Roll  rig  installed  in  IAR  6x9  ft  and  AFRL  7x10  ft  wind  tunnels 
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Fig.  3  Pitch  rig  and  its  installation  in  USAF  7x10  ft  wind  tunnel 


Fig.  5  65°  delta  wing  model 


6a  slow  dissipation  6b  fast  dissipation 

Fig.  6  Different  spiral  breakdown  appearances 
(A=65°,  a=26°,  <>=0o) 
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Fig.  9  Response  of  vortex  breakdown  location  to  perturbations  (a=36°) 


Fig.  10  Response  of  vortex  breakdown  location  to  perturbations  (a=34°) 


Fig.  1 1  Bifurcation  of  vortex  breakdown  locations 
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Fig.  14  Rolling  amplitude  effect  on  vortex  breakdown 
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Fig.  15  Measured  vortex  breakdown  location  during  ramp-and-hold  motion  in  roll 
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Fig.  15(Cont.)  Measured  vortex  breakdown  location  during  ramp-and-hold  motion  in  roll 


Fig.  16  Breakdown  movement  during  pitch-down  (a:30°— >20°) 


Fig.  17  Breakdown  movement  during  pitch-up  with  different  rates  (oc:20°— >50°) 
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18a.  0=24.43°  18b-  <*=24.40° 


Fig.  18  Breakdown  fluctuation  during  pitch-down  motion  (a:50°— >0°) 


£2=0.51,  a=37.5°,<|)=2.50 
Fig.  19  flow  visualization  about  coning  effect  on  vortex  breakdown 
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Fig.  20  Coning  motion  effect  on  vortex  breakdown  movement 
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21c  coning  induced  leading-edge  shape  2 Id  coning  induced  camber  surface 

Fig.  21  Motion  induced  kinematic  changes  in  wing  geometry 


22b  Q  =0.62,  a=25  °,  <j)=0o 


22c  £2  =-0.62,  a=25°,  <])=0° 


Fig.  22  Leading-edge  vortex  trajectory  under  non-planar  motion 
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Fig.  23  Vortex  breakdown  locations  over  flat  wing  and  cambered  wing  halves  (())  =0) 
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Fig.  24  Estimated  and  measured  breakdown  locations  for  rolled  65°  delta  wing 
(g=30°) 
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Question  by  Dr.  Greenwell:  How  was  the  vortex  core  size  determined? 

Answer:  The  vortex  core  diameter  was  determined  by  applying  image  analysis  software 
to  images  of  section  view  of  the  vortex  ring.  The  gray  level  is  proportional  to  the  number 
of  particles,  which  is  the  indication  of  vorticity  strength.  The  location  of  maximum  slope 
in  gray  level  is  approximately  the  edge  of  vortex  core.  Details  can  be  found  in  my  paper 
at  the  8th  ISFV. 

Comment.  Dr.  Greenwell:  See  my  later  Paper  5:  The  radius  of  the  smoke  ring  can  be 
approximated  as 

I'Miiok,'  —  ti  (dp/3)  sc] it [ Re  (Xp—  1 )] 

Therefore,  rsmoke  is  an  indicator  of  the  vortex  swirl  ratio,  <f),  and,  therefore,  is  proportional 
to  the  vortex  strength,  F. 

Question  by  Prof.  Slooff:  Has  the  prediction  method  been  used  for  other  configurations 
than  the  one(s)  tested  by  the  authors  and,  in  particular,  has  it  been  applied  to  wings  with  a 
rounded  leading  edge? 

Answer:  We  have  used  this  method  to  predict  delta  wings  with  different  sweep  back 
angles  in  static  conditions.  A  preliminary  comparison  has  also  been  conducted  for  a  70° 
delta  wing  in  dynamic  (pitch)  conditions.  So  far  comparisons  are  only  for  sharp  leading- 
edge  delta  wings.  But  in  principle,  this  method  can  be  applied  to  other  shapes  of  leading 
edges  for  delta  wings. 
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